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We discuss preliminaries of a calculation of the form factors for the semileptonic decays B —> nlv, 
B s Klv, and B K£ + £~. We simulate with NRQCD heavy and HISQ light valence quarks on 
the MILC 2 + 1 dynamical asqtad configurations. The form factors are calculated over a range of 

{SJ momentum transfer to allow determination of their shape and the extraction of \V„b | ■ Additionally, 

we are calculating ratios of these form factors to those for the unphysical decay B s — s- tj^. We are 
studying the possibility of combining these precisely determined ratios with future calculations 

f^i of B s —> rj s using HISQ fe-quarks to generate form factors with significantly reduced errors. 
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1. Motivation 

We are improving upon our previous B — > nlv calculation [1] in several ways, including the 
use of: Zj-quark smearing; HISQ light valence-quarks with random wall sources [2]; better scale- 
determination [3]; fitting advances (e.g. simultaneous fits to multiple separation times); and the 
2-expansion [4]. The calculation will also benefit from improved experimental data [5] which, 
when combined with lattice results, determines |V,,/,|. 

In parallel, we are studying the B s — > Klv decay. In combination with planned measure- 
ments [6], this will provide an additional exclusive determination of \V u b . Not yet studied on the 
lattice, this decay has a heavier spectator quark than B — > niv and should have reduced errors. 

We are also studying B — > K£ + £~, where the flavor-changing neutral current b — > s provides a 
probe of new physics (cf. Ref. [7]). There are existing [8] and promised [9] experimental results 
for this decay, but few unquenched lattice calculations [10]. 

Additionally, we are investigating the possibility of using the unphysical B s — >• r\ s decay to 
build ratios of form factors using NRQCD ^-quarks in which the leading sources of error largely 
cancel. This ratio could be combined with a future calculation of B s — > T] s using a HISQ ft-quark to 
yield form factors with greater precision, ie. 



f(B -> xlv) 



X f(B s ^Tl s )\ mSQb , (1.1) 

NRQCD b 



f(B S -»• T]s) 

analogous to the recent HPQCD work on B and B s decay constants [11]. 
2. Calculation 

The Standard Model (V — A)* 1 weak interaction responsible for the b — >• u transition results in 
hadronic matrix elements (XjV^lB^), parameterized via form factors 

R Y „ / Ml -Ml \ n y Mn -Ml 

(X\V*\B q ) =fi X (a 1 ) U+rf— L^fl* +f B Q « X {q 2 ) \ 2 * (2.1) 

\ HI H 



where = pg — p%- We recast these form factors in terms of lattice-convenient form factors, 



(X\V»\B q ) = J2M Bi 



PB q J,X, 2n , H f B q X ( 2 N 



(2.2) 



where p^= Px~ {Px • Pb^bJm^ . In the B^-meson rest frame, the form factors are simply related 
to the temporal and spatial components of the hadronic vector matrix elements, 



(X\V°\B q ) = JlAhJ^iq 2 ) 

(X\V k \B q ) = x [2MB~ qP k x f B /{q 2 ). (2.3) 



We calculate the components of the hadronic vector matrix elements and, from them, construct the 
form factors f1 qX (q 2 ) for the decays listed in Sec. 1. Ultimately, the form factors are related to 
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ensemble 


a [fm] 


m/(sea)/m s (sea) 


Aronf 


AW 


L 3 xN t 


T 


CI 


0.12 


0.005/0.05 


1200 


2 


24 3 x 64 


12-15 


C2 


0.12 


0.01/0.05 


1200 


2 


20 3 x 64 


12-15 


C3 


0.12 


0.02/0.05 


600 


2 


20 3 x 64 


12-15 


Fl 


0.09 


0.0062/0.031 


1200 


4 


28 3 x 96 


21-24 


F2 


0.09 


0.0124/0.031 


600 


4 


28 3 x 96 


21-24 



Table 1: Left to right: ensemble, lattice spacing, light and strange sea-quark mass, number of con- 
figurations, number of source times, volume, and separation between parent and daughter mesons. 



experimentally measured differential decay rates 1 



dT B * x Gl\V ub \ 



dq 2 1927T 3 M 3 



M^+Ml-q 2 ) -AMI M 2 X 



2 



3/2 



\fT(q 2 )\\ (2.4) 



where experimental and lattice results are combined to determine Iv^ . The Standard Model sup- 
pressed b — > s transition in B — > K£ + l~ opens the door for potentially discernible new physics 
contributions. The search for new physics in this decay requires the tensor form factor, related to 
the (kO) -component of the hadronic tensor matrix element 



/ K \T k0 \B) = 2MbPk ifi K (q 2 ). (2.5) 
x 1 17 M B +M K JT KH ' v ' 

2.1 Generating Correlator Data 

Ensemble averages are performed using the MILC 2+1 asqtad gauge configurations [12] listed 
in Table 1. The valence quarks in our simulation are NRQCD [13] ^-quarks, tuned in Ref. [11], and 
HISQ [14] light and strange quarks, whose propagators were generated in previous works [2, 15]. 
Working in the parent meson rest frame, a sequential propagator is built from NRQCD b and 
spectator HISQ quarks. The ft-quark smearing function (y' — y) is either a delta function or 
Gaussian, specified by indices a,/3 in Eqs. (2.6, 2.8), and is introduced by the replacement £ y — > 
Eyy 0(y' — y)- The spectator source includes a U(l) phase E, (x'). The daughter quark, with U(l) 
phase and momentum insertion at x, is tied to the sequential quark propagator, with £ x in Eqs. (2.6 
- 2.8) accomplished via random wall sources, ie. £ x — > £ x x / 1, (x)^(x'). 



= 73 E<*S,&y) (2-6) 

x,y 

c x (t ,f,p) = ^£ e !>(x - y) (^(f,y) *U*o>*)) ( 2 - 7 ) 

L x,y 

Cg qX (t ,t,T;p) = 1 £ e<^-*Hd> x (t + T,x) J(t,z) ^(t ,y)) (2.8) 



x,y,z 



'Eq. (2.4) neglects final state lepton masses. 
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<P(y' - y) : 



: £(x'), C(x) 



t + T, x' 




Figure 1 : Setup for three-point correlator data generation. 




1-1, s-s 1-1, 1-s, s-1, s-s 



f (0) 



(a) Fits to aMg' for various source-sink combinations. 

0.552 




9 10 



(c) flMj"' vs. the number of states N, with fmm/a = 2. 




(b) Fits to b 

0.545 

0.544 
0.543 
0.542 



1.(0) 



1-1, 1-s, s-1, s-s 
for various source-sink combinations. 




(d) aM^ versus '»„„/", with N = 10. 



Figure 2: From ensemble C2: (a) and (b) show the improvement from simultaneous fits to local 
("1") and smeared ("s") source-sink combinations; (c) and (d) display stability of fit results. 



In three-point correlator data the parent meson is created at time-slice to, the daughter meson is 
annihilated at to + T, and a flavor-changing current J(t,z) is inserted at intermediate times to <t < 
to + T, where to is chosen at random to reduce auto-correlations. This three -point correlator setup 
is depicted in Fig. 1 . Data are generated over the ranges of parent and daughter meson temporal 
separations listed in Table 1 . Prior to fitting, all data are shifted to a common to = 0. 
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Figure 3: Fit results for E % and M K are 
combined with simulated pion momentum 
to check the dispersion relation. A plot of 
(E\ —M^)/p 2 vs. (rip) 2 is consistent with 
the expected result of 1 + <??(ap) 2 . 
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Figure 4: For B — > ntv on ensemble C3, fit results for 
a(Vj.) cont (defined in Sec. 2.3) are shown for combina- 
tions of T used in simultaneous fits. Colored bands 
correspond to the "best-fit" combination T = 14, 15. 



2.2 Fitting Correlator Data 

Two-point correlator data for parent mesons are fit to the ansatz 



n=0 



r (-iyv 



Bq 



where b, 



a(n) 
'1 



9 / 



2M. 



(n) 



(2.9) 



Data are generated and analyzed for the B and B s mesons and all four combinations of local and 
Gaussian smeared sources and sinks. Fig. 2 shows the improvement observed from simultaneously 
fitting multiple source-sink smearing combinations and the stability of fit results with respect to 
changes in Af and t mm . Two-point correlator data for the daughter mesons are fit to 



Cx(t) 



N-l 

I 

n=0 



I I4 n) l 2 (-ir 



where 



IE 



(n) 



(2.10) 



We generate and analyze data for the n, K, and r\ s daughter mesons, each at momenta p G 1% jh x 
{(0,0,0), (1,0,0), (1, 1,0), (1, 1, 1)}. These fit results satisfy the dispersion relation as shown in 
Fig. 3. At each momentum, three-point correlator data are fit to 

C« qX (t, T) = bf n) A^f 4 n)t e -< ) ' e -4°(r-0 ) ( 2.ii) 

m,n=0 

where the thiee-point amplitude is related to the lattice matrix element by 



(2.12) 



We perform a simultaneous, Bayesian fit to the four local and smeared combinations of the parent 
two-point, the daughter two-point, and three-point correlator data sets for multiple values of T . The 
improvement from simultaneously fitting data for multiple T is shown in Fig. 4. 
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(b) f +fi (q 2 ) for B s ->K£v. 




20 21 

q 2 |GeV 2 ] 

(c) f +fi (q 2 ) and if T (q 2 ) for B -> 




19 20 21 

q 2 |GeV 2 | 

(d)/ +i0 (^ 2 ) forfij -> Tjj. 



Figure 5: Preliminary form factor results. 

2.3 Matching and Preliminary Results 

The lattice vector current (/ = "V^) is matched to the continuum at one-loop using massless 
HISQ lattice perturbation theory [1, 16] 



,(l),subv 



(2.13) 



where (-^ (1) ' sub ) = (f^') - a^io^^TO- Currents contributing through &(a s ,*QCD/M,a>/aM) are 



(0) 



and 



(i) 



1 

2M 



bypf-Vq 



(2.14) 



where q is the daughter quark in Fig. 1. For the lattice tensor current (/ = =^ v )> 

(Tlo)cont = (i + a sPT ){3r$) + <5# ) ' Bub > (2.15) 

where (■3^ ) ' mh ) = OfcCio^ 05 )- Heavy-quark symmetry of the NRQCD fr-quark allows 

the tensor current renormalization to be recast in terms of vector current quantities: JTq = 'f, 



kO 



(I) 



k ' 



■7^, and Cfo = "W- 
For the ensembles analyzed, preliminary results for form factors are shown in Fig. 5. The form 
factors f+fl(q 2 ) are calculated for all decay channels and ifj(q 2 ) is calculated for B — > K£ + £~. 
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3. Next Steps 

Once data generation and correlator fitting is complete, we will extract the physical values of 
the form factors. The kinematic dependence of the form factors over the full range of physical 
q can be written in a model-independent way via the z-expansion [4]. We plan to incorporate 
the chiral and continuum extrapolations in the z-expansion as in [2, 15]. The resultant modified z- 
expansion permits the use of data over the full range of q 2 , including data at daughter momenta for 
which chiral perturbation theory is expected to break down. We plan to cross-check these results 
against those obtained by separately performing the chiral and continuum extrapolation and then 
the z-expansion. 
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